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ABSTRACT 

Statistical weak lensing by large-scale structure - cosmic shear - is a promising cos- 
mological tool, which has motivated the design of several large upcoming surveys. 
Here, we present a measurement of cosmic shear using coadded Sloan Digital Sky 
Survey (SDSS) imaging in 168 square degrees of the equatorial region, with r < 23.5 
and i < 22.5, a source number density of 2.2 per arcmin^ and median redshift of 
Zmed — 0.52. Th ese coadds were g enerated using a new method described in the com- 
panion Paper I (|Huff et al.ll201l[ ) that was intended to minimise systematic errors in 
the lensing measurement due to coherent PSF anisotropics that are otherwise preva- 
lent in the SDSS imaging data. We present measurements of cosmic shear out to an- 
gular separations of 2 degrees, along with systematics tests that (combined with those 
from Paper I on the catalogue generation) demonstrate that our results are domi- 
nated by statistical rather than systematic erro rs. Assuming a cosmological model 
corresponding to WMAP7 (jKomatsu et al.ir201l[ ) and allowing only the amplitude of 
matter fluctuations as to vary, we find a best-fit value of as — ■GSG'to'l'^i (If); with- 
out systematic errors this would be as = 0.636to^°37 (1^)- Assuming a fiat ACDM 
model, the combined constraints with WMAP7 are as = 0J8Ato!^Qll{la)toml{2a) 
and flmh^ = 0.1303+Q;oQ4g(lcr)l:Q;Q09(2o-); the 2a error ranges are respectively 14 and 
17 per cent smaller than WMAP7 alone. Aside from the intrinsic value of such cosmo- 
logical constraints from the growth of structure, we identify some important lessons 
for upcoming surveys that may face similar issues when combining multi-epoch data 
to measure cosmic shear. 

Key words: cosmology: observations ~ gravitational lensing: weak - surveys. 



1 INTRODUCTION 

As a result of gravitational lensing, large scale inhomo- 
geneities in the matter density field produce small but sys- 
tematic fluctuations in the sizes, shapes, and fluxes of dis- 
tant objects that are coherent across large scales. This effect 
was first suggeste d as a tool for constrainin g the form of the 
metric in 1966 bv lKristian fc Sachsl (| 19661 ). In a more mod- 
ern context, the two-point statistics of lensing fluctuations 
allow the only truly direct measurement of the matter power 
spectrum and the growth of structure at late times, when 



dark energy has caused an accelerated expansi on of the uni- 
verse iess et ahll 19981 : iPerlmutt^ et al.lll999l ') and affected 
the growth of structure. Many studies have pointed out 
that high signal-to-noise cosmic shear measurements would 
be extraordinarily sensitive probes of cosmological param- 
eters (e.g., Hiutcrcr 1998; Benabed & van W acrbckc 2004j), 
which led to its being flagged as one of the most promis- 
ing probes of dark en ergy by the Dark Energy Task Force 
dAlbrecht et al.|[2006l ). Direct measurements of the growth 
of structure also offer the opportunity t o test alternative 
models of gravity (e.g.. iLaszlo et al]|201ll '). 



© 0000 RAS 



2 E. M. Huff et al 



C osmic shear m easurements were attempted as early as 
1967 (|Kr istian" 1967), but until the tu r n of the millennium 
dBaconet al. 2000: Kaiser et al.ll2000l : Ivan Waerbeke et all 
l2000l : IWittman et al.ll200Gl ). no astronomical survey had the 
statistical power to detect it. The difficulty of the measure- 
ment is a consequence of the near- homogeneity and isotropy 
of the universe. An order unity distortion to galaxy images 
requires an integrated line-of-sight matter over-density of: 



ds 



AnG cLl dLS 



(1) 



where ds, dL, and dLS are the angular diameter distances 
from the observer to the background source, from the ob- 
server to the lens, and from the lens to the background 
source, respectively. A fluctuation in the surface density AE 
leads to a shear distortion 7 ~ AE/Ecrit. 

Averaged over large (~ 100 Mpc) scales, typical line-of- 
sight matter fluctuations are only 10~'^Ecrit. The primary 
source of noise in the shear measurement, the random in- 
trinsic dispersion in galaxy shapes, is orders of magnitude 
larger; typically the shape noise results in a dispersion in 
the shear of a-y = 0.2. Worse, even in modern ground- 
based astronomical imaging surveys, the coherent distor- 
tions - or point-spread function (PSF) - induced by effects 
of the atmosphere, telescope optics, and detectors are typi- 
cally several times la rger than the cosmological signal (e.g., 
iHevmans et al.ll201ll and Paper I in this series). Estimat- 
ing t he distances to the backgroun d sources is both cru- 
cial jMa et al.l I2OO6I) and difficult (|Ma fc BernsteinI I2OO8I : 



iBernstein &: Huterej I2OIOI '): errors there will modulate the 
amplitude of the signal through Ecrit, biasing inference of 
the growth of structure. 

These obstacles define the observational problem. While 
the existence of cosmic shear has been established by 
the first studies to detect the effect, the full potential of 
cosmological lensing remains to be exploited. Few data 
sets capable of achieving the signal strength for a cos- 
mologically competitive measurement presently exist - the 



Canada-France-Hawaii Telescope Legacy Survey (CFHTLS; 


Hoekstra et al. 


l2006l:ISemboloni et al.ll2006l: iBeniamin et all 


20071: Fu et al. 


20081). the Cosmological Evolution Survey 



and the subset of the SDSS imaging studied here. How- 
ever, several large surveys are planned for the immediate 
and longer-term future that will substantially expand the 
amount of available data for cosmological weak lensing stud- 
ies. In the n ext few years, these include Hyper Suprime- 
Cam f HSC. iMivazaki et al.1 120061). Dark Energy S urvey 
(DE^I, iThe Dark Energv Survev Collaboration! |2005| ). the 
Kilo-Degree Survey (KID^J, and the Panoramic Survey 
Telescope and Ra pid Response System (Pan-STARRqf|, 
iKaiser et al.|[201ol ). Further in the future, there are even 
more ambitious prograr ns such as the Large Synoptic Sur- 
vey Telescope (LSSlH ILSST Science Collaborations et af] 



[iooi), Euclicfl, and the Wide-Field Infrared Survey Tele- 
scope (WFIRSllfl). 

For this work, we have combined several methods dis- 
cussed in the literature as viable techniques for measuring 
cosmic s hear while remov ing common systematic errors. In 
Paper I (|Huff et al.ll201ll ). we began with the PSF model 
generated by the Sloan Digital Sky Survey (SDSS) pipeline 
over ~ 250 deg^ that had been imaged many times, and em- 
pl oyed a rounding kernel m ethod similar to that proposed 
in IBernstein fc Jarvid (|2002l ). The result, after appropriate 
masking of problematic regions, was 168 square degrees of 
deep coadded imaging with a well controlled, homogeneous 
PSF and sufficient galaxy surface density to measure a cos- 
mic shear signal. The usable area in r band was only 140 
square degrees because of a PSF model error problem on the 
camcol 2 charge-coupled device (CCD), which is suspected 
to be an amplifier non-linearity problem. 

In this work, we use the catalogue from Paper I to 
produce a cosmic shear measurement that is dominated by 
statistical errors. Section |3] enumerates the primary sources 
of systematic error when measuring cosmic shear using our 
catalogue (the properties of which are summarized briefly in 
Sec.[5|, and describes our approaches to constraining each of 
them. In Section |4l we outline our correlation function esti- 
mator and several transformations of it that are used for sys- 
tematics tests. Our methods for estimating covariance ma- 
trices for our observable quantities (both due to statistical 
and systematic errors) are described in Sec. [5] Finally, sec- 
tion [S] presents the constraining power of this measurement 
alone for a fiducial cosmology, and in combination with the 
7-year Wilkinson Mi crowave Anisotropy Probe (WMAP7, 
iKomatsu et al.l I2OIII ) parameter constraints to produce a 
posterior probability distribution over Q,mh? , ^bh? , o"8, ns, 
and w. We show that in addition to its value as an inde- 
pendent measurement of the late-time matter power spec- 
trum, this measurement provides some additional constrain- 
ing power over WMAP7 within the context of ACDM. We 
conclude with some lessons for the future in Sec. [7] 

While this work was underway, we learned of a parallel 
effort by Lin et al. (2011). These two efforts use different 
methods of coaddition, different shape measurement codes, 
different sets of cuts for the selection of input images and 
galaxies, and analyze their final results in different ways; 
what they have in common is their use of SDSS data (not 
necessarily the same sets of input imaging) and their use 
of the SDSS Photo pipeline for the initial reduction of the 
single epoch data and the final reduction of the coadded 
data (however, they use different versions of Photo). Us- 
ing these different methods, both groups have extracted the 
cosmic shear signal and its cosmological interpretations. We 
have coordinated submission with them but have not con- 
sulted their results prior to this, so these two analysis efforts 
are independent, representing versions of two independent 
pipelines. 



^ https : //www . darkenergy survey . org/ 
^ http : //www. astro-wise . org/projects/KIDS/ 
http : //pan-Starrs . if a.hawaii . edu/public/ 
* http://www.lsst.org/lsst 



^ http : //sci . esa. Int/science-e/www/area/index . cf m?f areaid 
=102 

^ http://wfirst.gsfc.nasa.gov/ 
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2 CATALOGUES 



Paper I l|Huff et al.ll201ll ) describes a coadd imaging dataset, 
optimised for cosmic shear measurement, constructed from 
single-epoch SDSS images in the Stripe 82 equatorial region, 
with right ascension (RA) —50° <RA< +45°, and decli- 
nation — 1.25° <Dec< +1.25°. In that work, we apply an 
adaptive rounding kernel to the single-epoch images to null 
the effects of PSF anisotropy and match to a single homo- 
geneous PSF model for the entire region, and show that in 
the resulting shear catalogues, the amplitude of the galaxy 
shape correlations due to PSF anisotropy at angular sepa- 
rations greater than 1 arcminute is negligible compared to 
the expected cosmic shear statistical errors. 

The final shape catalogue described in that work con- 
sists of 1 067 031 r-band and 1 251 285 i-band shape measure- 
ments with characteristic limiting magnitudes of r < 23.5 
and i < 22.5, over effective areas of 140 and 168 square 
degrees, respectively. 



amplitudes of T^-p^f^psf, ^int, (7eint), the uncertainty in the 
theoretically-predicted ^-y,_B arising from redshift errors, and 
the uncertainty in the shear calibration (via the responsivity 
TV) contribute less than 20 per cent to the statistical errors 
in (ee). 

Our approach to handling of systematic error is as fol- 
lows: we attempt to reduce each systematic to a term that 
can be robustly and believably estimated from real data (ei- 
ther the data here or in other, related work), and we then ex- 
plicitly correct for it. These corrections naturally have some 
uncertainty associated with them, which we use to derive a 
systematic error component to the covariance matrix. The 
exception to the rule given here is if there is a systematic 
error for which there is no clear path to estimating its mag- 
nitude, then we do not attempt any correction, and simply 
marginalize over it by include an associated uncertainty in 
the covariance matrix. 



3.1 Cosmic shear 



3 MODEL FOR THE LENSING AND 
SYSTEMATIC ERROR SIGNALS 

We model the observed galaxy shape field as the sum 
of a cosmic shear component, an independent systemat- 
ics field produced by anisotropies in the effective PSF 
Bpsf, and a systematics field produced by the intrinsic 
spatial cor relations of galaxy sh apes eint (intrinsic align- 
ments; e.g.. lHirata fc Seliakll2004l ). We allow for a shear cal- 
ibration factor that depe nds on the shear responsivity TZ 
l|Bernstein fc Jarvisll200 of the ensemble of galaxy surface 
brightness profiles to the underlying gravitationally-induced 
shear 7. Typically 7?. ~ 1 — e^^s, however we consider it to 
be a more general factor that also includes any biases due 
to effects such as uncorrected PSF dilution, noise-related bi- 
ases, or selection biases. We assume that the galaxy shape 
response to PSF anisotropies TZpsi is not a priori known, but 
rather suffers from a similar set of 'calibration' uncertainties 
as the response of the ensemble of galaxy images to gravita- 
tional lensing shear. Thus we define our model for the two 
ellipticity components e — (61,62) as 



e = TZj + 7?.psf fipsf + Bin 



(2) 



We assume that the two-point statistics of the under- 
lying (cosmological) shear field (77) consist entirely of E- 
modes, e-y^s (which is a good enough approximation given 
the s ize of our errors: [Crittenden et aI.|[2002 : ISchneider et al.l 
l2002t ). and are statistically independent of the PSF when 
averaged over large regions. We also assume that the PSF 
and the intrinsic alignments are independent - but not that 
the lensing shear and intrinsic alignments are independent 
l|Hirata fc SeliakI |2004| ). The two-point correlation of the 
galaxy shapes contains terms resulting from gravitational 
lensing and from systematic errors: 



(ee) = ATZ^^^^E + 7?.psfCpsf + ^im + (7eint) 



(3) 



Here, ^paf is the auto-correlation of the PSF ellipticity field. 
Errors in the determination of the galaxy redshift distribu- 
tion will enter as a bias in the predicted £,-y,E- 

Our goal is to carry out a statistics-limited measure- 
ment of ^■y,E- This will entail showing that the combined 



Foreground anisotropies in the matter distribution along the 
line of sight to a galaxy will generically distort the galaxy 
image. For weak lensing, the leading order lensing contribu- 
tion to galaxy shapes can be thought of as arising from a lin- 
ear transformation of the image coordinates Aa;true = a^obs , 
where 



1 + K + 71 
72 



72 

1 + K — 71 



(4) 



The convergence k causes magnification, whereas the 
shear components 71 and 72 map circles to ellipses. The 
shear is related to the projected line-of-sight matter distri- 
bution, weighted by the lensing efficiency: 



(71,72) = d~ 



W (x, XO {dl - dl, 2d.dy) 5 (xftO dx.(5) 



Here we integrate along the comoving line-of-sight distance 
X (where Xi is the distance to the source), and the matter 
over-density 5 = {p — ]}) /p. The window function in a fiat 
universe is 



W(X, X^) = ^fi™J^O (1 + 2)X' f - - - 

2 Vx Xi 



(6) 



The two-point correlation function of the shear can be 
calculated by identifying pairs of source galaxies, and defin- 
ing shear components (7*, 71) for each one to be the shear 
in the coordinate system defined by the vector connecting 
them, and in the 7r/4 rotated system. This two-point corre- 
lation function can be expressed as a linear transformation 
of the matter power spectrum Ps averaged over the line of 
sight to the sheared galaxies: 



and 



{7t7t) ± (7x7x) 



2n 



(7) 



\2d% 



dx 



d{x) 



' , / / '\ d{x' -xV^ 
dx (X j 



dix' 



(8) 
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where the last expression makes use of Limber's approx- 
imation and d(x) is the distance function, i.e. x in a 
flat universe, K~^^^ siuK^^^x in ^ closed universe, and 
{—K)~^^'^ sinh(— Ar)^''^X in ^n open universe. In the expres- 
sion in brackets, n(x') represents the source distribution as 
a function of line-of-sight distance (normalised to integrate 
to 1). This statistic (P^) is sensitive both to the distribution 
of matter 5 and to the background cosmology, via both the 
explicit n„i dependence and the distance-redshift relations. 



3.2 Intrinsic alignments 

Many studies have discussed intrinsic alignments of galaxy 
shapes due to effects such as angular momentum alignments 
or tidal torque due t o the large-scale densi t y field (for pi- 
oneering studies, se e ICroft fc MetzlejbOOOl: iHeavens et al.l 



20001: ICatelanetal 



l200ll : ICrittenden et al.ll200ll : |jing||20ol 



Hopkins et al.ll2005f ). While these effects can generate co- 
herent intrinsic alignment 2-point functions, [Hir ata & S eliakI 
l|2004 ) pointed out that the large-scale tidal fields that can 
cause intrinsic alignments are sourced by the same large- 
scale structure that is responsible for producing a cosmic 
shear signal. Thus, in this model, the intrinsic alignments 
do not just have a nonzero auto-correlation, they also have 
a significant anti-correlation with the lensing shear which 
can persist to very large transverse scales and line-of-sight 
separations. If left uncorrected, this coherent alignment of 
intrinsic galaxy shapes suppresses the lensing signal, since 
the response of the intrinsic shape to an applied tidal field 
has the opposite sign from the response of the galaxy image 
to a shear with the same magnitude and direction. We gener- 
ally refer to the intrinsic alignment auto-correlation as the 
"77" contamination and its correlation with gravitational 
lensing as the "G7" contamination. This can be compared 
to the pure gravitational lensing auto-correlation ("GG"). 

To address the uncertainty related to intrinsic 
alignments, we rely on empirical measurements that 
constrain the degree to which they might affect our 

using SPSS imag i ng a nd 



we 

the degree 
measurement. Several studies 



spectroscopic dat a (e.g., iMandelbaum et al.l l2006a 



Hirata et al. '2007"; 'Okumura e t al.l l2009l : I Joachimi et al.l 
2011 : Mandolbaum et al. 2Qllb|) have demonstrated the 
existence of intrinsic align ments of ga l axy s hapes on 
cosmological distance scales. iHirata et al.l (|2007l ) used the 
luminosity and colour-dependence of intrinsic alignments 
for several SDSS galaxy samples to estimate the contamina- 
tion of the cosmic shear signal due intrinsic alignments for 
lensing surveys as a function of their depth. These estimates 
were a function of the assumptions that were made, for 
example about evolution with redshift. The "central" model 
given in that paper leads to a fractional contamination of 



-0.08 



(9) 



for a limiting magnitude of m_R,iim = 23.5, which is close 
to the limiting magni tude of ou r sample. Subseq uent work 
l|Joachimi et all 1201 ll : [Mandelbaum et all boilbl ) provided 
more information about redshift evolution; primarily those 
results were in broad agreement with the previous ones, and 
were sufficient to rule out both the "op t imisti c" and the 
"very pessimistic" models in lHirata et al.l (|2007f ). 

We thus adopt the "central' model, and apply the cor- 



rection given in Eq. ((9| to our theory predictions for the Ct 
due to cosmic shear, multiplying the predicted cosmic shear 
power spectrum by 0.92 before transforming into the statis- 
tics that are used for the actual cosmological constrainti|3|- 
We also assume this correction has a conservative system- 
atic uncertainty of 50 per cent, which amounts to an overall 
4 per cent uncertainty in the theory prediction (see Sec. [5] 
for a quantitative description of how we incorporate this and 
other systematic uncertainties into the covariance matrix). 

Since the G7 correlation is first order in the intrinsic 
alignment amplitude, while the 77 power is second order, we 
expect the first to be the dominant systematic. In principle, 
the G7 effect could be smaller than 77 if the correct align- 
men t model is quadratic in the tidal field rather than lin- 
ear (|Hirata fc Seliakll2004l ). However, in the aforementioned 
cases in which intrinsic alignment signals are detected at 
high significance (i.e. for bright ellipticals) the linear model 
for in trinsic alignments appears to be valid (|Blazek et al.l 
I2OIII) . Therefore we attempt no correction for 77. 



3.3 Shear calibration 

Another source of systematic error for weak lensing mea- 
surements is uncertainty in the shear calibration factor. The 
galaxy ellipticity (e+,ex) observed after isotropizing the 
PSF need not have unit response to shear: in general, aver- 
aged over a population of sheared galaxies, we should have 



((e+,ex)) = 7^(7+,7x), 



(10) 



where TZ is the shear responsivity. It depends on both the 
shape measurement method and the galaxy popul ation (e.g. 
iMassev et a"l]|2007bl : lBernstein|[201ol : IZhandboill ) . 

For this work, we used the re-Gaussianization method 
jHirata fc Seliakll2003l ). which is based on second moments 
from fits to elliptical Gaussians, and has been previously 
applied to SDSS single- epoch imaging ([Mandelbaum et al.l 
I2OO5I : iReves et ahlbOllf ). For this class of methods, in the 
absence of selection biases and weighting of the galaxies, 
perfectly homologous isophotes, and no noise , there is an 
analytic expectation l|Bernstein fc Jarvisll200 j ): 



7^ = 2(l 



(11) 



where Crms is the root-mean-square ellipticity per component 
(+ or x). 

The calibration errors for re-Gaussianization and other 
adaptive-we ighting methods are well-studied in the liter- 
ature (e.g., iHirata et all l2004bl : IMandelbaum et al.l |2005| . 
l2011al : lReves et allboill ). They arise from all of the devia- 
tions from the assumptions of Eq. . Higher-order depar- 
tures from non-Gaussianity in the galaxy light profile cause 
errors in the PSF dilution correction. Errors in the mea- 
surement of the PSF model will cause a similar error in the 
dilution correction. The resolution factor of an individual 
galaxy depends on its ellipticity, so any resolution cut on 
the galaxy sample will introduce a shear bias in the galaxy 



^ While the intrinsic alignm ent contam ination is in principle 
scale-dependent, the plots in Hirata et al.l |2007i) suggest that this 
scale dependence is in fact quite weak for the scales used for our 
analysis, so we ignore it here. 
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selection function. Due to the non-linearity of the shear in- 
ference procedure, noise in the galaxy images causes a bias 
in the shears (rather than just making them noisier). The es- 
timation of the shear responsivity, or even of erms, is another 
potential source of error, as the response of the galaxies to 
the shear depends on the true, intrinsic shapes, rather than 
the gravitationally sheared, smeared (by the PSF), noisy 
ones that we observe. 

Past approaches to this problem have included detailed 
accounting for these effects one by one. In this paper, we 
instead use detailed simulations of the image processing 
and shape measurement pipelines, including real galaxy im- 
ages, to estimate both the shear calibration and the red- 
shift distribution of our catalogue. The advantage is that 
this includes all of the above effects and avoids uncertain- 
ties associated with analytic estimates of errors. The Shera 
(SHEar Reconvolution An alysis) simulation packagafi has 
been previously described (jMandelbaum et al.l l2011al ) and 
applied to single-e poch SDSS data for galaxy-galaxy lensing 
l|Reves et al.ll201ll ). but this is its first application to cosmic 
shear data. 

To simulate our images, we require a fair, flux-limited 
sample of any galaxies that could plausibly be resolved in our 
coadd imaging, including high-resolution images with realis- 
tic morphologic^. For this purpose we use a sample of 56 662 
galaxy images dra wn from the COSMOlo g ical Evolution 
Survey (COSMOS: iKoekemoer et al. I l2007l : IScoville et al.l 
l2007a lbl) imaging catalogues. The deep Hubble Space Tele- 



scope (HST) Advanced Camera for Surveys/ Wide Field 
Camera (ACS/WFC) imaging in F814W ("broad /") in this 
1.6 deg^ field is an ideal source of a fairly-selected galaxy 
sample with high resolution, deep image^ij. These images 
consist of two samples - a "bright" sample of 26 116 galax- 
ies in the magnitude range I < 22.5, and a "faint" sample 
consisting of the 22.5 < I < 23.5 galaxies . The charge trans- 
fer ineffi ciency-corrected ( Massev et al.l |2010|) and multi- 



drizzled (jKoekemoer et al 



2OO2I : [Rhodes et all [ioOTl . to a 



pixel scale of 0.03") galaxy postage stamp images have been 
selected to avoid CCD edges and diff'raction spikes from 
bright stars, and have been cleaned of any other nearby 
galaxies, so they contain only single galaxy images without 
image defects. The b right sample is used for gro und-based 
image simulations in iMandelbaum et al.l (|2011al ) ; the faint 
sample is selected and processed in an identical wa-yf"! Each 
postage stamp is assigned a weight to account for the rel- 
ative likelihoods of generating postage stamps passing all 
cuts (avoidance of CCD edges and bright stars) for galax- 
ies of different sizes in the COSMOS field; this weight is 
calculated empirically, by comparing the size distribution of 
galaxies with postage stamps to the size distribution of a 
purely fiux-limited sample of galaxies. 



* http : //www. astro .prlncet on. edu/~rmandelb/shera/ shera. html 
^ Simple models with analytic radial profiles and elliptical 
isophotes are not adequate to measure all sources of systematic 
error such as under-fitting biases or those due to non-elliptical 
isophotes {Bernstein 201Ci). 

Admittedly there may be some sampling variance that affects 
the morphological galaxy mix. 

We thank Alexie Leauthaud for kindly providing these pro- 
cessed images. 



Each of these postage-stamp images has several proper- 
ties associated with it that are of interes t for this analysis . 
The COSMOS photometric catalogues (jllbert et al.1 120091 ) 
contain HST _F814H^ magnitudes as well as photometric 
redshifts and Subaru r — i colours based on PSF-matched 
aperture magnitudes. 

In order to simulate our observations, we first select a 
coadd 'run' consisting of five adjacent frames in the scan 
direction at random from the list of completed runs. We 
draw 1250 galaxies (exactly 250 per frame) at random from 
the list of COSMOS postage stamps according to the weights 
described above, up- weighting the probability of drawing the 
faint galaxies by a factor of 1.106 to account for the fact that 
we have sampled the faint population at a lower rate than 
the bright one in constructing the image sample. 

Once a list of postage-stamp images is selected, we as- 
sign r- and i-band magnitudes by re-scaling each image; each 
galaxy image is inserted into the coadded imaging with the 
flux it would have been observed to have in SDSS before 
the addition of pixel noise. The i-band is chosen to be 0.03 
magnitudes brighter than the COSMOS F81AW (I) band 
MAG_AUTO values; this small offset is based on empirical com- 
parison with SDSS magnitudes for brighter galaxies, to ac- 
c ount for slight differences in the FSIAW and i passbands 
(Ma ndelbaum et al.ll2011a[ ). The r-band is chosen so as to 
match the Subaru PSF-matched aperture colours for each 
object. Each postage stamp is assigned a random, uniformly- 
sampled position in the coadd run, with the postage stamps 
distributed equally among the frames. 

We use the SHERA code to pseudo-deconvolve the HST 
point-spread function, apply (if necessary; see below) a shear 
to each galaxy, reconvolve each image with the known coadd 
point-spread function, renormalise the flux appropriately, 
and resample from the COSMOS pixel scale to the coadd 
pixel scale before adding that posta ge stamp to t he coadd 
image. This procedure, sugg est ed bvlKaiseil ||2000|) and im- 
plemented to high precision in IMandelbaum et al.l (|2011al ). 
can be used to simulate ground-based images with a shear 
appropriately applied, despite the space-based PSF in the 
original COSMOS images, and with a user-defined PSF. 

The normal coadd masking algorithm is then applied, 
and shear catalogues are generated as in Paper I by run- 
ning the SDSS object detection and measurement pipeline. 
Photo-frames, followed by the shape measurement code 
described in Sec. 13.31 The output catalogues are matched 
against the known input object positions, and a simulation 
catalogue of the matches is created. We employ these simu- 
lations below to determine the shear calibration and as an 
independent validation of our inferred redshift distribution. 

For each suite of simulation realisations, we use the 
same random seed (i.e., we select the same galaxies from 
our catalogue and place them at identical locations in the 
coadded image) but with different applied shears per com- 
ponent ranging from —0.05 to -f 0.05. We measure the mean 
weighted shape of the detected simulation galaxies produced 
by our pipeline, and fit a line to the results. Since the same 
galaxies are used without rotation, only the slope and not 
the intercept is meaningful. The shear response in each com- 
ponent for each applied shear is shown in Fig. [T] The re- 
sponsivities in the two components are consistent, which is 
expected on oversampled data with a rounded PSF. (The 
unequal size of the error bars reflects the number of runs 
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Shear calibration simulations 
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Figure 1. The response of the mean ellipticities (ei) and (e2) 
to applied shear, as determined in the SHERA-based simulations. 
Poisson error bars are shown. The additive offset to the response 
curve is not shown in the fit; these simulations do not accurately 
measure an additive shear bias. 



that we were able to process by the time the shear calibra- 
tion solution was frozen.) The total number of galaxies in 
the final simulated catalogues was 130 063. The response ap- 
pears to be linear for small applied shears. Based on these 
results, we adopt a shear responsivity for this galaxy popu- 
lation of 1.776 ± 0.043. 



3.4 Redshift distribution 

The explicit dependence of the shear signal in Eqs. ((5} 
and ([8} on the distribution of lensed galaxy redshifts, com- 
bined with the practical impossibility of acquiring a spectro- 
scopic redshift for the millions of faint galaxies statistically 
necessary for a cosmic shear measurement, can be a trouble- 
some source of bias and systematic uncertainty for cosmic 
shear measurements. 

An error in the estimated redshift distribution leads 
to an incorrect prediction for the amplitude of the shear 
signal at a given cosmology. This is similar in principle to 
the bias arising in the amplitude of the galaxy-galaxy lens- 
ing signal due t o pho tometric redshift biases explored in 
iNakaiima et al.l l|201 j ): uncorrected, standard photometric 
redshift estimation techniques can lead to biases in the pre- 
dicted lensing signal at the ~ 10 per cent level. For cosmic 
shear measurements, an imperfect estimate of the redshift 
distribution leads to biases in ag and fim that are compara- 



ble in amplitude to t he errors in the estimated mean of the 
redshift distribution llvan Waerbeke eral]|2006D . 

As a fiducial reference, the redshift distribution of the 
single-epoch SDSS ima ging catalogue is es tablished to ap- 
proximately 1 per cent jSheldon et al.|[201ll ): for deeper sur- 
veys over a smaller area, this becomes a more difhcult prob- 
lem, as the spectroscopic calibration samples available for 
inferring the redshift distribution are limited in their red- 
shift coverage and widely dispersed across the sky. We em- 
plo y a colour-matching technique similar to that employed 
by Sheldon et al.l l|20lj ^: in what follows, we describe the 
technique, our estimate of its uncertainty, and several cross- 
checks on the results. 



3.4..! Fiducial redshift distribution 

The source red shift distribution used in our analysis is de - 
rived following IHma et al I ll200i) and ICunha et al.l l|2009f ). 
and is similar in spirit to lSheldon et alT ( 2011 ): the princi- 
ple is that, for two galaxy samples that span broadly similar 
ranges in redshift, colour, and limiting magnitude, matched 
colour samples correspond to matched redshift distributions. 

Our spectroscopic calibration sample is composed of 
12 360 galaxie s, from th e union of the VIMOS VLT Deep 
Survey (Le F evre et al.i r2005. VVDS) 22h field , the DEEP2 
Gala xy Redshift Survey ( Davis et al.ll2003l : iMadgwick et al.l 
l2003f ). and portions of th e PRism MUlti-object Survey 
(PRIMUS: [Coil et al.|[20Tll . Cool et al. 2011 in p rey.) We 
follow the procedures outlined in INaka iima et al. (2011) for 
selecting good quality spectroscopic redshifts, and avoiding 
duplicate galaxies in samples that overlap (such as DEEP2 
and PRIMUS). Each of these samples has a redshift distri- 
bution that is likely to differ substantially from the redshift 
distribution of our lensing catalogue: the DEEP2 catalogue 
in the fields we use at 23'' 30" and 02'' 30" is heavily colour- 
selected (in non-SDSS bands) towards objects at z > 0.7; the 
PRIMUS catalogue includes several fields, some of which are 
selected from imaging with a shallower limiting magnitude; 
and the VVDS catalogue is selected in the I band (/ < 22.5) 
with a relatively high redshift failure rate that exhibits some 
colour dependence. 

We assign a redshift from a galaxy in the union calibra- 
tion sample to the closest galaxy in the lensing catalogue 
within 3 arcsec, finding 12 360 matches. To generate a rep- 
resentative training sample of galaxies from the lens cata- 
logue, we draw 4 x 10^ galaxies with replacement from the 



full area (not just in these regions), with sampling proba- 
bility proportional to the mean of the weights assigned in 
the r and i bands to that galaxy for the correlation analy- 
sis (Eg. I17|) . Note that this procedure does not incorporate 
those galaxies i n the exclud e d cam col 2 region. 

We use the iLima et al.l (|2008l ) cod43 to solve for a set 
of weights over the calibration sample, such that the re- 
weighted 5-dimensional magnitude distributions of the cal- 
ibration sample match those of the representative random 
subset of the lensing catalogue. The histogram of the cali- 
bration sample redshifts reweighted in this manner is shown 
as a solid line in Fig. [21 The inferred mean redshift is 0.51; in 



http : //kobayashi .physics . Isa.umich . edu/~ ccunha/nearest/ 
© 0000 RAS, MNRAS 000, 000-000 



Cosmic shear in SDSS 7 



ir Matching Technique, <z> = 0,51 
COSMOS simulotions, <z> - 0.50 
SDSS single-epoch, <z> - 0.42 




(13) 



To propagate these errors into the covariance matrix for 
(,E, we first fit a smooth function of the form 



riz (z) cc z e ^ ' ' 



(14) 



to the nearest neighbour weighting-derived redshift dis- 
tribution shown in Figure [2] the best fit parameters are 
a = 0.5548, zo = 0.7456, and b = 2.5374. We perturb 
this smooth distribution by adding a random number drawn 
from a normal distribution with mean (zj ) (normalised to 
the weighted number of cahbration galaxies in that bin) and 
standard deviation aj at the location of the jth redshift bin. 
We then renormalise the perturbed distribution to unity, and 
compute the predicted cosmic shear signal. The covariance 
matrix of 402 realisations of this procedure is added to the 
statistical covariance matrix. 



Figure 2. The redshift distribution inferred from matching the 
colours of the spectroscopic calibration sample to those of the 
lensing catalogue (solid black line, Sec. 13.4. Il l shown alongside 
the noisier redshift distribution inferred from the shear calibration 
simulations (dashed red line, Sec. l3.4."3ll . The bes t-fit distribution 
for th e single-epoch SDSS lensing catalogue from Naka iima et al.l 
ll201lh is shown for reference as the blue dot-dashed line. 

contrast to the redshift distribution for single-epoch imag- 
ing, there is a non-negligible fraction of the galaxy sample 
above z > 0.7. We use the solid curve based on the colour- 
matching techniques to calculate the shear covariance ma- 
trix, and to predict the shear correlation function for any 
given cosmology. 



3.4-2 Uncertainty 

We expect that the primary source of error in the red- 
shift distribution as estimated from the combined calibra- 
tion sample is sample variance, resulting from the finite vol- 
ume of the calibration s ample. To est i mate its magnitude, 
we use the public code of lMoster et al.l (|201ll ) for estimating 
the cosmic variance of number counts in small fields. 

Our redshift binning scheme has 19 bins between < 
z < 1.5. For a collect i on of disparate calibration fields, we 
use the iMoster et al.l (|201ll ) code to produce a fractional 
error in the number counts cTgg^ij- for the jth redshift bin in 
the i field (where fields are distinguished by their coverage 
area) in bins of stellar mass. 

The redshift sampling rate of each distinct survey in 
the calibration sample differs, and so the balance of contri- 
butions to the final redshift distribution will change as well. 
To account for this, we sum over every calibration field's 
contribution to the reweighted redshift distribution in the j 
bin to estimate an absolute (not relative) overall error: 

= i'^SS,i,jncH,i,jf (12) 
i 

where the effective number of galaxies contributed in the j 
bin by the i survey is just the sum over the nearest-neighbour 
derived weights assigned to calibration sample galaxies k in 
that field i and bin j: 



3.4.3 Other tests 

As an independent check on the redshift distribution, we 
also use the shear calibration simulations (Sec. I3.3|l to con- 
strain the redshift distribution of our sources. The COSMOS 
photometric redshifts, inferred as they are from many more 
imaging bands (typically with deeper imaging) than for the 
SDSS data discuss ed here, are very accurate. For example, 
lllbert et al.l l|2009h finds a photo- 2: scatter of o"z/(l -I- z) ~ 
0.01 for a galaxy sa mple with the flux lim it of the SDSS 
coadds. In contrast, iNakai ima et al.l ()201lh found that in 
the SDSS single-epoch imaging, the scatter defined in the 
same way was ~ 0.1 despite the brighter flux limit of the 
single-epoch imaging (due in part to the more limited num- 
ber of bands, but primarily to the far lower signal-to-noise 
ratio). If we treat the COSMOS photometric redshifts as 
we would spectroscopic data, then the redshift distribution 
of COSMOS galaxies that pass successfully into the shear 
catalogue is the same as that of our source catalogue - as- 
suming, of course, that the COSMOS fleld is representative 
of the whole of Stripe 82. It is not, of course; large-scale 
structure in the COSMOS fleld (which c an be signiflcant, a s 
COSMOS covers only 1.7 square degrees: iKovac et al.ll201oh 
can bias a determination of the redshift distribution in this 
manner. The n{z) inferred from the COSMOS-based sim- 
ulations is also shown in Fig. O and agrees extremely well 
with the fiducial n{z) derived from colour-matching. 

A final (but obviously not independent) s anity check is 
to compare to the COSMOS Mock Catalogue ijjouvel et al.l 
;2009), which is being used extensively to plan future dark 
energy programmes, using the cuts reff > 0.47 arcsec, limit- 
ing magnitudes r < 23.5, and i < 22.5 (see Paper I, where 
we argue that these most closely mimic the cuts in our data). 
This test predicts (z) = 0.51, identical to that obtained via 
the re-weighting procedure. Given the crudeness of the pro- 
cedure for comparing the results, this is an excellent val- 
idation of the COSMOS Mock Catalogue as a forecasting 
tool. 

3.5 Stellar contamination 

Stellar contamination of the galaxy catalogue reduces the 
apparent shear by diluting the signal with round objects that 
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are not sheared by gravitational lensing. Because the image 
simulations described in Sec. l3.3l onlv included galaxies, the 
resulting shear responsivities do not include signal dilution 
due to accidental inclusion of stars in the galaxy sample. In 
Paper I, we estimated the stellar contamination by compar- 
ison with the DEEP2 target selection photometry (which is 
deeper and was acquired at the Canada- France-Hawaii Tele- 
scope under much better seeing conditions than typical for 
SDSS), and found a contamination fraction of 0.017. We also 
argued that the mean stellar density in the stripe must be 
larger than in the high-latitude DEEP2 fields, by a factor 
as large as 2.8. We therefore conservatively take the stellar 
contamination fraction /star to be 



Mean ellipticity: r band 



/star = 0.017(1.9 ± 0.9) = 0.032 ± 0.015. 



(15) 



The resulting suppression of the cosmic shear signal is 
treated in much the same way as for intrinsic alignments: we 
reduce the theory signal by a factor of (1 — 0.032)^ = 0.936, 
and add a contribution to the covariance of 0.030 times the 
theory signal. 



3.6 Additive systematics 

Among the most worrying systematics in the early detec- 
tions of cosmic shear was additive power. This comes from 
any non-cosmological source of fluctuations in shapes such 
as PSF anisotropy that add to the ellipticity correlation 
function of the galaxies. Such power was clearly detected 
in Paper I in the form of systematic variation of both star 
and galaxy ei as a function of declination. The sense of 
the effect - a negative contribution to ei (in r band we 
havcEl (ei) = -0.0018 and (ea) = -1-0.0004, while in i band 
(ei) = -0.0022 and (ea) = -0.0002) - is suggestive of mask- 
ing bias, in which the selection of a galaxy depends on its 
orientation, with galaxies aligned in the along-scan direction 
(ei < 0) being favoured, and with no effect on 62 (consis- 
tent with zero mean over the whole survey). The reason 
for this particular sign is seen in Figure 2 of Paper I; as 
shown, bad columns along the scan direction tend to be re- 
peated at the same location in multiple images, resulting 
in significant (non-isotropic) masks with that directional- 
ity. Direct evidence for masking bias comes from the change 
in mean ellipticity due to increased masking: when we re- 
moved from the coadded image pixels that were observed in 
fewer than 7 input runs and reran Photo-Frames, the (ei) 
signal became worse: —0.0051 in r band and —0.0044 in i 
band, whereas (62) was essentially unchanged. This increase 
is difficult to explain in terms of spurious PSF effects, so we 
conclude that our galaxy catalogue likely contains a mixture 
of masking bias as well as possible additive systematics from 
PSF ellipticity in the coadded image. 

The mean ei signal as a function of declination is shown 
in Fig.Oin bins of width 0.05 degrees. We take this as a tem- 
plate for mask-related selection biases (combined with any 
systematic uncorrected PSF variation as a function of dec- 
lination, which in west-to-east drift-scan observations is a 
highly plausible type of position dependence). Before com- 



The la Poisson uncertainty in these numbers is 0.0005 (0.0004) 
per component in r (i) band. 




-0.5 0.5 
Declination (deg) 

IVIean ellipticity: 1 band 



1.5 




-0.5 0.5 1 1.5 
Declination (deg) 

Figure 3. The mean ellipticity (ei) as a function of declination 
in the r and i bands. This signal was removed from the galaxy 
catalogue prior to computing the final correlation function. The 
r band data between declination —0.8° and —0.4° were rejected 
due to the known problems with camcol 2. The error bars are 
Poisson errors only. 

puting the correlation function, we subtracted this mean sig- 
nal from the galaxy ellipticity catalogue. 

One danger in this procedure to remove spurious (ei) 
is that some real power could be removed - that is, even 
in the absence of any systematic error, some of the actual 
galaxy shape correlation function signal could be suppressed 
since the method determines the mean ei of the real galax- 
ies and by subtracting it introduces a slight artificial anti- 
correlation. The best way to guard against this is with simu- 
lations. Using the Monte Carlo simulation tool of Sec. 15.1.21 
we generated simulated realisations of our ellipticity cata- 
logue and either implemented the (ei) projection or not. 
The difference in the correlation functions is a measure of 
how much power was removed. The result is shown in Fig. HJ 
and shows that the loss of real power is insignificant com- 
pared to our error bars. 



3.6.1 PSF anisotropy 

Convolution with an elliptical PSF will induce a spuri- 
ous ellipticity in observed galaxy surface-brightness profiles. 
While the effective PSF for these coadds is a circular double 
Gaussian to quite high precision, the tests in Paper I indi- 
cate a low level of residual anisotropy that we must consider 
here. 
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Figure 5. The star-galaxy ellipticity correlation functions. Shown are the rr, ri (i.e. star r X galaxy i), ir, and ii correlation functions, 
reduced to 10 bins. The solid points, which are offset to slightly lower (9-values for clarity, are the H — h correlation functions, and the 
dashed points are the X X functions. All error bars arc Poisson only. 
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Figure 4. The loss of actual power due to ei projection. Us- 
ing 36 realizations from the Monte Carlo simulation, we find the 
difference in post-projection ellipticity correlation function ^(6) 
and original 5(9). These are shown as the solid points (5+-|-) and 
dashed points (5x x ) in the figure, re-binned to 10 bins in angular 
separation 6. The dashed lines at top and bottom arc the ±lcr 
statistical error bars of our measurement. The reduction of actual 
power is detectable by combining many simulations, but is very 
small compared to the error bars on the measurement. 



Possible sources of this issue include: (i) inaccuracies 
in the single-epoch PSF model used to determine the ker- 
nel to achieve the desired PSF; (ii) colour-dependence of 
the PSF that means the single-epoch PSF model from the 
stars is not exactly the PSF for the galaxies; or (iii) the 
fact that we determine the rounding kernel on a fixed grid, 
so that smaller-scale variations in PSF anisotropy might re- 
main uncorrected. All of these must be present at some level, 
although the last two cannot be the full solution: (ii) does 
not explain the residual stellar ellipticit-vF*l. and (iii) does 
not explain why there is structure in the declination direc- 
tion on the scale of an entire CCD (0.23 degrees). 

For a galajcy and a PSF that are both well- 
approximated by a Gaussian, the PSF- correct ion given 
above produces a measured ellipticity of: 



obs PSF 1 — i?2 PSF 
e = /-Cpsfe — 6 



R2 



(16) 



see e.g. iBernstein fc Jarvisl (|2002l ). The weighted (by the 

We have searched for a, g — i dependence in the stellar elliptic- 
ities in the coadded image. We only found effects at the ~ 0.002 
level, and while they are statistically significant, we have not es- 
tablished whether they correspond to true colour dependence ver- 
sus e.g. variation of stellar colour distributions along the stripe. 
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Contamination implied by star-galaxy correlation 
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Figure 6. The implied contamination to the galaxy ellipticity 
correlation function if the star-galaxy correlation function is used 
as a measure of the additive PSF power. The solid points are 
the H — h correlation functions, and the dashed points are the X X 
functions. All error bars are propagated from the Poisson errors 
assuming correlation coefficient -1-1 (a better assumption than in- 
dependent errors, but likely an overestimate). The dotted curves 
show the 1(T errors in each radial bin from the Monte Carlo sim- 
ulations (see Sec. 15. 1.21 1 which include both Poisson and cosmic 
variance uncertainties. Note also that the shapes and normalisa- 
tions of the H — |- and X X signals are nearly identical. 

same weights used for the correlation function; see Eq. I17|) 
average of the PSF anisotropy response defined in Eq. ^ 
over the sample of galaxies considered in this work is TZ^sf ~ 
0.86 (r band) or 0.95 {i band); in what follows we take a 
value of 0.9. 

A nonzero star-galaxy correlation function resulting 
from systematic PSF anisotropy (as estimated in Paper I) in- 
dicates the presence of a spurious contribution to the shear- 
shear correlation function with amplitude « Q.Q^sg. We will 
not determine this response to high enough accuracy to sub- 
tract the effect with small residual error: doing so would not 
require just a simulation, but a simulation that knows the 
correct radial profile of the PSF errors0 In our case, the 
star-galaxy correlation function is detectable but below the 
errors on the galaxy-galaxy ellipticity auto-correlation (al- 
though not by very much), so a highly accurate correction 
is unnecessary. 

We constrain the PSF anisotropy contribution by com- 
puting the star-galaxy correlation function. This was done 
in Paper I, but some of the star-galaxy signal is due to the 
systematic variation of PSF ellipticity with declination and 
is removed by the subtraction procedure above. The star- 
galaxy ellipticity correlation function with the corrected cat- 
alogue is shown in Fig. [5] The implied contamination to the 
galajcy ellipticity correlation function, appropriately averag- 
ing the bands and applying the factor of 7?.psf = 0.9, is shown 
in Fig. [6l 

These measured star-galaxy correlations can be used 



This might be an option in future space-based surveys if 
the type of error can be traced to the source of ellipticity 
(astigmatismxdefocus, coma, or jitter). In either space or ground- 
based data, one could imagine doing cross-correlations of higher- 
order shapelet modes l||lcfregier 20o3) to extract the particular 
form of the errors. None of these options are pursued here. 



to construct a reasonable systematics covariance matrix for 
this systematic. We take the amplitude of the diagonal ele- 
ments of the PSF systematic covariance to be equal to the 
amplitude of the measured contamination. We also assume 
that the off-diagonal terms are fully-correlated between bins, 
which is equivalent to fixing the scaling of this systematic 
with radius, and saying that only the overall amplitude of 
the systematic is uncertain. 

Since there are a number of uncertainties in this pro- 
cedure, we do not apply any correction for these additive 
PSF systematics as we do for ones that are previously dis- 
cussed, such as intrinsic alignments or stellar contamina- 
tion. Instead, we simply include a term in the systematics 
covariance matrix to account for it. We also will present a 
worst-case scenario for the impact of this term on cosmo- 
logical constraints; in Sec. [5] we will show what happens to 
the cosmology constraints if we assume that the systematic 
error is +2a from its mean, i.e. 40 per cent of the statistical 
errors. This should be taken as a worst-case scenario for this 
particular systematic. 



4 ANALYSIS TOOLS 

4.1 Ellipticity correlation function 

We compute the ellipticity correlation functions defined 
in Eq. (O on scales from 1-120 arcminutes. For the cos- 
mological analysis, we start by computing the correlation 
function in 100 bins logarithmically spaced in separation 
6 to avoid bin width artifacts. For the cosmological pa- 
rameter constraints, we project these onto the Complete 
Orthogonal Sets of E -/B-mode Integrals (COSEBI) basis 
(|Schneider et al.ll2010l ) to avoid the instabilities of inverting 
a large covariance matrix estimated via Monte Carlo simu- 
lations (we will describe our implementation of COSEBIs in 
Sec. 14. 3|) . However, for display purposes, it is more conve- 
nient to reduce the 9 resolution to only 10 bins so that the 
real trends are more visually apparent. 



4.1.1 Weighting 

The correlation functions used here are weighted by the in- 
verse variance of the ellipticities, where the "variance" in- 
cludes shape noise. Specifically, we define a weight for a 
galaxy 



1 



(t| + 0.372 ' 



(17) 



where is the ellipticity uncertainty per component de- 
fin ed by our shape m easurement pipeline. As demonstrated 
bv lReves et al.l (|201ll ). these may be significantly underesti- 
mated in certain circumstances; however, this will only make 
our estimator slightly sub-optimal, so we do not attempt to 
correct for it. The value of 0.37 for the root-mean-square 
(RMS) intrinsic el lipticity dispersion per component comes 
from the results of lReves et al.l (|201ll ). for r < 22, and there- 
fore we are implicitly extrapo lating it to fainter magnitudes. 
Given that iLeauthaud et aU ((2007 ) found a constant RMS 
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Figure 7. The ellipticity correlation functions in the rr, ri, ii and ww (combined) band combinations. The sohd points denote the H — h 
and the dashed points denote the X X components of the correlation function. The points have been slightly displaced horizontally for 
clarity. The Monte Carlo errors are shown. 



ellipticity to far fainter magnitudes in the COSMOS data, 
we consider this extrapolation justifietCfl. 



4-1-2 Direct pair- count code 

A direct pair-count correlation function code was used for 
the cosmological analysis. It is slow (~ 3 hours for 2 x 10® 
galaxies on a modern laptop) but robust and well-adapted to 
the Stripe 82 survey geometry. The code sorts the galaxies 
in order of increasing right ascension a\ the galaxies are 
assigned to the range —60° < a < +60° to avoid unphysical 
edge effects near q = 0. It then loops over all pairs with 



N ote that we do not use the actual value of RMS ellipticity 
from [Leauthaud et al.l ||2007| ) - only the trend with magn itude 
— because, as demonstrated bv Mandelbau m et al.l l l2011al ), the 
RMS ellipticity value inXeauthaud et al. (20071) is not valid for 
our adaptively-defined moments, which use an elliptical weight 
function matched to the galaxy light profile. 



|ai — 02] < Smax. The usual ellipticity correlation functions 
can be computed, e.g. 



(18) 



where the sum is over pairs with separations in the rele- 
vant 6 bin, and the ellipticity components are rotated to 
the line connecting the galaxies. The direct pair-count code 
works on a fiat sky, i.e. equatorial coordinates (a, 5) are 
approximated as Cartesian coordinates. This is appropri- 
ate in the range considered, \5\ < 1.274°, where the maxi- 
mum distance distortions are i 

5max = 2.5 X 10"*. The direct 
pair-count code is applicable to either auto-correlations of 
galaxy shapes measured in a single filter (rr, ii) or cross- 
correlations between filters or between distinct populations 
of objects (ri and all of the star-galaxy correlations). 

Simple post-processing allows one to compute the 5+ 
and 5_ correlation functions, defined by 



C+(e) = C++(e) + exx(e) 



(19) 
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and 

^_(0) = C++W-?xx(e). 

4-1.3 Combining bands 



(20) 



Finally, the different band correlation functions rr, ri, and 
ii must be combined according to some weighting scheme: 



(21) 



where the label "ww" indicates that the bands were com- 
bined. The relative weights were chosen according to the 
fraction of measured shapes in r- and i-bands, i.e. lo, 



J r : 



Wri = 2frfi, and wu — where the weights are fr — 0.4603 
and fi = 0.5397. 

The final ellipticity correlation functions (with the 6 
resolution reduced to 10 bins) are shown in Fig. [71 



4.2 Tests of the correlation function 

We implement several null tests on the correlation function 
to search for remaining systematic errors. 

The first test, shown in Fig. [S] constructs the difference 
between the cross-correlation function of r and i band galaxy 
ellipticities versus the rr and ii auto-correlations. The dif- 
ferences in the two types of correlation functions are small 
compared to the statistical uncertainty in the signal. This is 
consistent with our expectations, as the true cosmic shear 
signal should be independent of the filters in which galaxy 
shapes are measured. 

The second test, shown in Fig. (9] compares the (band 
averaged or ww) correlation function computed using galaxy 
pairs separated in the cross-scan (north-south) direction ver- 
sus pairs separated in the along-scan (east-west) direction. 
This difference should be zero if the signal we measure is 
due to lensing in a statistically isotropic universe. The er- 
ror bars shown are Poisson errors, so they may be slight 
underestimates at the larger scales, where cosmic variance 
becomes important. Visual inspection shows no obvious off- 
set from zero, but the error bars are larger for this test than 
in Fig. [8] because the null test includes no cancellation of 
galaxy shape noise. 



4.3 E/B-mode decomposition 

As a final check for systematics, we decompose the 2-point 
correlation function (2PCF) into E- and _B-modes, where, 
to leading order, gravitational lensing only creates _B-modes. 
The B-modes c an arise from the li mited validity of th e Born 
approximation (I Jain et al. 20001: Hilbert et al.ll2009l ). red- 
shift source clusterin g (ISchneider et aL 2002h. and lensing 
(mag nification) bias (| Schmidt et al. 20091 : iKrause fc Hiratal 
I2OIOI ). however the amplitude of B-modes from these sources 
should be undetectable with our data. At our level of sig- 
nificance, a B-mode detection would indicate remaining sys- 
tematics, e.g. due to spurious power from an incomplete PSF 
correction. 

Formerly used methods to decompose i5-and B-modes, 
such as the aperture mass dispersion 
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Figure 8. The difference between the galaxy ellipticity cross- 
correlations (ri) and the auto-correlations [rr + ii) /2, with error 
bars determined from the Monte Carlo simulations. The upper 
panel shows the H — h correlations and the lower panel shows the 
X X correlations. The dashed line is the Icr statistical error bar 
on the actual signal. 



with the filter functions T± as derived in ISchneider et al.l 
(2002), or the shear E-mode correlation fu nction, suffer from 
_E/B-mode mixing (jKilbinger et al.ll2006l ). i.e. B-modes af- 
fect the i?-mode signal and vice versa. These statistics can 
be obtained from the measured 2PCF, for an exact E/B- 
mode decomposition, however they require information on 
scales outside the interval [Smin;&max] for which the 2PCF 

has been measured. 

The ring statistic s ('Schneide r fc Kilbingeij l2007l : 
lEifler et al.1 l2010l; IFu fcK ilbingcr 20ldi) and more recently 
the COSEBIs ("Sc hneider et a l. 2010) perform an EB-mode 
decomposition using a 2PCF measured over a finite angular 
range. COSEBIs and ring statistics can be expressed as 
integrals over the 2PCF as 



EB 
and 

i?EB(6l) = 



d6_ 
T 



(23) 



dff 
29 



-[i+[e')z+{e\e)±i.{e')z-{e',e)]. (24) 



For the ring statist i cs, we use the filter functions Z± spec- 
ified in lEifier et al] (|2010D . The de rivation of the COSEB I 
filter functions T±n is outlined in 'Schneider et al.1 (|2010t ). 
where the authors provide linear and logarithmic filter func- 
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Null test, [N/S]-[E/W] pairs: rr++ 



Null test, [N/S1-[E/W] pairs: ri++ 
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Figure 9. The null test of the correlation functions measureti using galaxy pairs whose separation vector is within 45° of the north-south 
direction, minus that measured using galaxy pairs whose separation vector is within 45° of the east-west direction. The error bars shown 
are the Poisson errors only. The dashed curve shows the lu error bars of the actual signal (all colour combinations and separation vectors 
averaged). The 6 panels show the three colour combinations (rr, ri, and ii) and the 2 components (H — h or X x). 



tions indicating whether the separation of the roots of the 
filter function is distributed linearly or logarithmically in 9. 
Note that whereas the ring statistics are a function of an- 
gular scale, the COSEBIs are calculated over the total an- 
gular range of the 2PCF, condensing the information from 
the 2PCF naturally into a set of discrete modes. The lin- 
ear T-functions can be expressed conveniently as Legendre 
polynomials, however compresses the cosmological in- 
formation into significantly fewer modes; we therefore choose 
the logarithmic COSEBIs as our second-order shear statis- 
tic in the likelihood analysis in Sec. |6l The COSEBI filter 
functions are displayed graphically in Fig. IIUI 

Figure [11] shows three different _E/i3-mode statistics de- 
rived from our measured shear-shear correlation function, 
i.e. the COSEBIs, the ring statistics, and the aperture mass 
dispersion. The error bars are obtained from the square root 
of the corresponding covariances' diagonal elements (statis- 
tics only). Note that the COSEBIs data points are signifi- 
cantly correlated. Slightly smaller is the correlation for the 
aperture mass dispersion, and the ring statistics' data points 
have the smallest correlation. 

From the COSEBIs, we find a reduced for the E- 
modes to be consistent with zero of 6.395, versus 1.096 for 
the B- modes (5 degrees of freedom each). The latter is con- 
sistent with purely statistical fluctuations. 



5 COVARIANCE ESTIMATION 

5.1 Ellipticity correlation function covariance 
matrix 

The covariance matrix of the ellipticity correlation function 
estimated via Eq. (|2ip was computed in several ways. The 
preferred method for our analysis is a Monte Carlo method 
(Sec. I5.1.2|l but we compare that covariance matrix with an 
estimate of the Poisson errors (Sec. IS.l.ip as a consistency 
check. 

5.1.1 Potsson method 

The direct pair-count correlation function code can compute 
the Poisson error bars, i.e. the error bars neglecting the cor- 
relations in ei+ej+ between different pairs. This estimate of 
the error bar is 

Equivalently, this is the variance in the correlation function 
that one would estimate if one randomly re-oriented all of 
the galaxies. The Poisson method is simple, however it is not 
fully appropriate for ri cross-correlations (since the same in- 
trinsic shape noise is recovered twice for pairs that appear 
in both ri and ir cross-correlations). Moreover, at scales of 
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Figure 11. The measured COSEBIs, ring statistics, and aperture mass dispersion from the combined cosmic shear signal. The error 
bars equal the square root of the corresponding covariances' diagonal elements (statistics only). Note that the COSEBIs data points are 
significantly correlated. Slightly smaller is the correlation for the aperture mass dispersion, and the ring statistics' data points have the 
smallest correlation. 
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Figure 10. The COSEBI filter functions T„+ (upper panel) and 
T„_ (lower panel) for the first 5 modes. 

tens of arcminutes and greater there is an additional contri- 
bution because the cosmic shear itself is correlated between 
pairs. Therefore the Poisson error bars should be used only 
as a visual guide: they would underestimate the true uncer- 
tainties if used in a cosmological parameter analysis. 

5.1.2 Monte Carlo method 

We used a Monte Carlo method to compute the covariance 
matrix of ^++(6) and ^xx(^)- The method is part theoret- 
ical and part empirical: it is based on a theoretical shear 



power spectrum, but randomizes the real galaxies to cor- 
rectly treat the noise properties of the survey. The advan- 
tages of the Monte Carlo method - as implemented here 
- are that spatially variable noise, intrinsic shape noise in- 
cluding correlations between the r and i band, and the sur- 
vey window function are correctly represented. The principal 
disadvantages are that the cosmic shear field is treated as 
Gaussian and a p articular cosmology must be assumed (see 
lEifler et al.ll2009l . for alternative approaches). However, so 
long as this cosmology is not too far from the correct one 
(an assumption that can itself be tested!), the Monte Carlo 
approach is likely to yield the best covariance matrix. 

The Monte Carlo approach begins with the generation 
of a suite of 459 realizations of a cosmic shear field in har- 
monic space according to a theoretical spectrum. For our 
analysis, the theoretical spectrum was that from the WMAP 
7-year (Larson ct al. 2011) cosmological parameter set (fiat 
ACDM; Q.bh^ = 0.02258; n^/i^ = 0.1334; ris = 0.963; 
Ho = 71.0kms~^ Mpc~^; and gg — 0.801), and th e shear 
power spectr um code used in [Albrecht et al.l (|2009l ) . itself 
based on the Eiscnstoin & Hu' (^1998') transfer function and 
the . Smith ct al. (2003) nonlinear mapping. The redshift dis- 
tribution discussed in section 13.4.11 based on a calibration 
sample from DEEP2, VVDS, and PRIMUS, was used as the 
input to the shear power spectrum calculation. 

From this power spectrum we generate a sample set 
of Gaussian i5-mode shear harmonic space coefficients a^. 
The full power spectrum is used at I ^ 1500; a smooth 
cutoff is applied from 1500 < / < 2000 and no power at 
/ ^ 2000 is included. This is appropriate for a covariance 
matrix since the power at smaller scales is shot noise dom- 
inated and cannot be recovered. (The i5-mode power spec- 
trum is CxMO = 3.6 X 10~^^, as compared to a shot noise 
of "{irit/n ~ 1.8 X 10~^.) No -B-mode shear is included. 
The particle - mesh spherical harmonic transform code of 
iHirata et all (|2004al ') with a 6144 x 3072 grid (L' = 6144) 
and a 400-node interpolation kernel {K — 10) was used to 
transform these coefficients into shear components (71,72) 
at the position of each galaxy j 

The use of a full-sky approach for the Monte Carlo realisations 
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A synthetic ellipticity catalogue was then generated as 
foUows. For each galaxy, we generated a random position 
angle offset -i/ij G [Oi 'L) '^^'^ rotated the ellipticity in both 
r and i bands by ^fjj o We then added the synthetic shear 
weighted by the shear responsivity to the randomised ellip- 
ticity to generate a synthetic ellipticity: 



Error bar comparison: MC/Poisson 



syn 2ii/jT true 

e/ =e ^^e, 



+ 1.737(ni 



(26) 



The 1.73 prefactor was estimated from Eq. which we 

expected to be good enough for use in the Monte Carlo 
analysis, so that the Monte Carlos could be run in parallel 
with the shear calibration simulations. The latter gave a final 
result of 1.78 ± 0.04, which is not significantly different. 

The direct pair-count correlation function code, in all 
versions (rr, ri, and ii) was run on each of the 459 Monte 
Carlo realisations, before combining the different correla- 
tions to get the weighted value via Eq. ((2T|. 

The Monte Carlo and Poisson error bars are compared 
in Fig. [T21 The correlation coefficients of the correlation 
functions in different bins are plotted graphically in Fig. 1131 

From each Monte Carlo correlation function we com- 
pute the COSEBls via Eq. ((23} and use their covariance 
matrix in our subsequent likelihood analysis. In order to 
test whether our covariance has converged, meaning that 
the number of realisations is sufficient to not alter cosmo- 
logical constraints, we perform 3 likelihood analyses in erg 
vs. Q.m space varying the numbers of realisations from which 
we compute the covariance matrix (see Sec. [6] for detailed 
methodology; for now we are just establishing convergence 
of the covariance matrix). In Fig. [14] we show the 68 and 
95 per cent likelihood contours, i.e. the contours enclose the 
corresponding fraction of the posterior probability (within 
the ranges of the parameters shown). We see that the con- 
tours hardly change when going from 300 to 400 realizations 
and show no change at all when going from 400 to 459 reali- 
sations, hence the 459 Monte Carlo realizations are sufficient 
for our likelihood analysis. 



5.2 Systematic contributions to the covariance 
matrix 

The following additional contributions are added to the 
Monte Carlo covariance matrix (and if appropriate the the- 
ory result) described in Sec. 15.1.21 

(i) The intrinsic alignment error was included following 
Sec. 13.21 the theory shear correlation function was reduced 
by a factor of 0.92, and an uncertainty of 4 per cent of the 
theory was added to the covariance matrix, i.e. we add an 
intrinsic alignment contribution 



Cov[^i, ^j] (intrinsic alignment) = 0.04 



2e(th)f(th) 



(27) 
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Figure 12. The ratio of error bars obtained by the Monte Carlo 
method to those obtained by the Poisson method, for 10 angular 
bins. The four curves show either rr or ii band correlation func- 
tions, and either the H — h or X X component. Note the rise in the 
error bars at large values of the angular separation, due to mode 
sampling variance. 
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Figure 13. The matrix of correlation coefRcicnts for the com- 
bined {ww) correlation functions in the 10 angular bins for which 
the correlation function is plotted in the companion figures. The 
bin number ranges from 0—9 for £^j^^(9) and from 10-19 for 
?x X (S); all diagonal components are by definition equal to unity. 
Based on 459 Monte Carlo realisations. 



where the theory curve (th) is obtained at the fiducial 



was not necessary for the SDSS Stripe 82 project, but was the 
simplest choice given legacy codes available to us. 

To simplify bookkeeping, the actual implementation was that 
a sequence of 10^ random numbers was generated, and a galaxy 
was assigned one of these numbers based on its coordinates in a 
fine grid with 0.36 arcsec cells in (a, 5). 



WMAP7 point. This covariance matrix includes perfect cor- 
relation between radial bins, implying that we treat this sys- 
tematic as being an effect with a fixed scaling with separa- 
tion, so the only degree of freedom is its amplitude. 

(ii) The stellar contamination was included following 
Sec. 13.51 the theory shear correlation function was reduced 
by a factor of 0.936, and an uncertainty of 3 per cent of the 
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Figure 14. Convergence test of the trg vs. f!m parameter con- 
straints as a function of the number of Monte Carlo reaUzations 
used to compute the covariance. The plot shows the 68 and 95 per 
cent likelihood contours (however, the lower 95 per cent contours 
are not visible). The covariance includes statistical errors only. 

theory was added to the covariance matrix, i.e. we add a 
stellar contamination contribution 

Cov[^i,,Cj] (stellar contamination) = 0.03^4'""'^]*', (28) 

where the theory curve (th) is obtained at the fiducial 
WMAP7 point. 

(iii) The implied error from the redshift distribution un- 
certainty is derived from 402 realisations of the sampling 
variance simulations as described in Sec. 13.4.2] We construct 
the covariance matrix of the predicted £-mode COSEBIs. 

(iv) The shear calibration uncertainty was conservatively 
estimated in Sec. 13.31 to be ±2.4 per cent, or equivalently 
4.8 per cent in second-order statistics. We thus add another 
term to the covariance matrix, 

Cov[Ci, CiKshear calibration) = 0.048^Cf '"^Cf"'- (29) 

(v) In Sec. 13.61 we described a procedure for including 
uncertainty due to additive PSF contamination. According 
to this procedure, the relevant systematics covariance matrix 
is related to the amplitude of the measured contamination 
signal: 

Cov[^i, 5j](PSF contamination) — 0.9^^sg,j^sg,i, (30) 

again assuming a fixed scaling with radius for this system- 
atic uncertainty. Since all entries scale together, we do not 
spuriously "average down" our estimate of the systematic 
error by combining many bins. 

The final data vector and its covariance matrix (includ- 
ing all the statistical and systematic components) are given 
in Tables IXTI and IX2l Note that given our procedure of ap- 
plying the systematic corrections to the theory, the data 
vector is the observed one without any such corrections for 
the stellar contamination and intrinsic alignments contami- 
nation. With this in hand, we can estimate the significance 



of the E- and B-mode signals described in section IT3] The 
probability that the COSEBI iJ-mode signal that we ob- 
serve is due to random chance given the null hypothesis (no 
cosmic shear) is 6.0 x 10~^. The probability of measuring 
our _B-mode signal due to random chance given the null hy- 
pothesis of zero B modes is .36, evidence that there is no 
significant B-mode power. 



6 COSMOLOGICAL CONSTRAINTS 

Having described the measured cosmic shear two-point 
statistics, and shown that the systematic bias in this mea- 
surement is small compared with the statistical constraints, 
we now turn to the cosmological interpretation. We work in 
the context of the flat A CDM parametrisatio n, taking where 
necessary the WMAP7 (|Komatsu et al.ll201lD constraints for 
our flducial parameter values. 

6.1 The prediction code: modeling second-order 
shear statistics 

To produce a cosmological interpretation of our measured 
cosmic shear signal from our model framework, we require 
a method to convert a vector of cosmological parameters 
into a prediction of the observed cosmic shear signal. Due to 
projection effects, we expect that a significant fraction of the 
observed cosmic shear signal is produced by the clustering of 
matter on nonlinear scales, so a suitably accurate prediction 
algorithm must ultimately rely on numerical simulations of 
structure formation. 

The prediction code used in our likeli hood analysi s is a 
modified version of the code described in lEifleil (|201ll ). We 
combine Haloflt (jSmith et al.ll2003l ). an analytic approach 
to modeli ng nonlinear structure , with the Coyote Universe 
Emulator (iLawrence et al.ll2010l ). which interpolates the re- 
sults of a large suite of high-resolution cosmological simu- 
lations over a limited parameter space, to obtain the den- 
sity power spectrum. The derivation is a two-step process: 
First, we calculate the linear power spectrum from an ini- 
tial power law spectrum Ps(k) oc fc"' e mploying the dewig- 
gled transfer function of lEisenstein fc H u (1998). The non- 
linear evolution of the density held is incorporated using 
Halofit. In order to simulate wC DM models we follow the 
scheme implemented in ICOSMO (|Refregier et al.ll201ll ). in- 
terpolating between flat and open cosmological models to 
mimic Quintessence cosmologies (see[Schrabback ot al. 201^ 
for more details). In a second step, we match the Halofit 
power spectrum to the Coyote Universe emulator (version 
1.1) power spectrum, which emulates Pg over the range 
0.002 fc 3.4/i/Mpc within ^ 2 ^ 1 to an accuracy of 
1 per cent. Wherever possible, the matched power spectrum 
exactly corresponds to the Coyote Universe Emulator; of 
course this is limited by the cosmological parameter space of 
the Emulator and its limited range in k and z. However, even 
outside the range of the Emulator, we rescale the Halofit 
power spectrum with a scale factor Pg (Coyote) /Pg (Halofit) 
calculated at the closest point in parameter space (cosmo- 
logical parameters, fc, and z) where the Emulator gives re- 
sults. Outside the range of the Emulator, the accuracy of this 
"Hybrid" density power spectrum is of course worse than 1 
per cent, however it should be a significant improvement 
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Figure 15. The 68 and 95 per cent likelihood contours of the 
combined data vector including a full treatment of systematics 
when using the Halofit prediction code (dashed) and when using 
the Coyote Universe-calibrated prediction code (solid). The red 
lines correspond to the best-fitting value of erg for a given • 
The dot indicates the WMAP7 best-fitting values. 

over a density power spectrum from Halofit only. From the 
so-derived density power spectrum we calculate tlie shear 
power spectrum via Eq. ^ and the shear-shear correlation 
function via Eq. (0. As a final step, we transform these 
predicted correlation functions to the COSEBI basis as de- 
scribed above in Sec. 14.31 

For our final results in the (Sim, as) likelihood analysis, 
we used both prediction codes; the results are compared in 
Fig. 1151 where they are seen to agree to much better than 
la. We therefore conclude that uncertainty in the theory 
predictions is sub-dominant to the other sources of system- 
atic error, and to the statistical error. 

6.2 Constructing the input data vector 

For our primary science results, we use the measured 5 
COSEBI modes (see Fig. 1111 left panel). As a first step we 
want to determine the number of COSEBI modes that need 
to be included in our likelihood analysis. In Fig.[TS]we show 
a likelihood analysis in the ag-ilm parameter space varying 
the number of modes in the data vector. We find that there 
is hardly a change in the likelihood contours when going 
from 4 to 5 modes; we therefore conclude that 5 modes is 
a sufficient number to capture the cosmological information 

encoded in our data set. 

As shown in lEifler et al.l (|2008l ) , the information con- 
tent of the aperture mass dispersion can be greatly improved 
when including 1 data point of the shear-shear correlation 
function ^+ into the data vector; here we adopt this concept 
for the COSEBIs. The basic idea is that the data point of 
the correlation function is sensitive to scales of the power 
spectrum to which the COSEBIs are insensitive. We incor- 
porate only a single data point of the correlation function 
as this is sufficient to capture the bulk of the additional in- 
formation while simultaneously minimising possible B-mode 
contamination. 



Figure 16. Convergence test of the ag vs. Qm parameter con- 
straints as a function of number of COSEBI modes in the data 
vector. The plot shows the likelihood contours enclosing 68 and 
95 per cent of the posterior distribution. (The lower bounding 
curve for the 95 per cent contours is not visible on the plot.) The 
covariance contains statistical errors only. The dot indicates the 
WMAP7 best-fitting values. 

In order to determine the optimal scale of the data point 
that is to be included, we consider 10 bins of f+ ranging from 
1.3 to 97.5 arcmin and perform 10 likelihood analyses for a 
combined data vector consisting of 5 COSEBI modes and 
one additional data point of We quantify the information 
content through the so-called q figure of merit (g-FoM) 

1 ~ \/|Q! , where 

Qij = y d^7rp(7r| d) (ttj - 7r-)(7rj - TTj) , (31) 

TT — {flm, as) is the parameter vector, p(7r|d) is the posterior 
likelihood at this parameter point, and tt' denotes the fidu- 
cial parameter values. If the likelihood in parameter space 
(i.e. the posterior probability) is Gaussian, the q-FoM cor- 
responds to the more common Fisher matrix based figure of 
merit / — 1/\\/F\. The Fisher matrix F can be interpreted 
as the expectation value of the inverse parameter covariance 
evaluated at the maximum likelihood estimate parameter 
set, which in our ansatz corresponds to the fiducial param- 
eters. Mathematically we can express this equivalence as 

f^^^^/\c7\^^\=q. (32) 
V I* 

Since the assumption of a Gaussian posterior is clearly vio- 
lated in the as-^m parameter space, we perform a full like- 
lihood analysis and calculate q to quantify the size of the 
likelihood. Note that smaller g-FoM is "better." 

We varied the angular scale (in arcmin) of the added 
(,+{0) data point, and found a minimal q-FoM a,t 6 — 37.8 
arcmin. We will use this scale for the additional ^+ data 
point henceforth. Note that this analysis uses a simulated 
input data vector in order to avoid biases from designing a 
statistical test based on the observed data. The constraints 
coming from the various possible data vectors - the COSE- 
BIs, the COSEBIs supplemented with a single ^+ point, and 
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COSEBIs 




Figure 17. The likelihood contours of the combined data vec- 
tor (solid), the shear-shear correlation function (dashed), and the 
COSEBIs (dotted) data vector to illustrate how much informa- 
tion is gained when including the additional data point. Note that 
the COSEBIs' lower 95 per cent contour is outside the considered 
region. The dot indicates the WMAP7 best-fitting values. 
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Figure 18. The correlation matrix of the COSEBI modes 1—5 
("E1...E5" in the figure) and ^+(38'). The left panel shows only 
the statistical (Monte Carlo) errors, and the right panel includes 
the systematics as well. 



the full shear correlation function - are compared in Fig. 1171 
They are not identical, which is expected since they weight 
the data in different ways, but are consistent with each other. 

The COSEBI modes are highly correlated with each 
other, and they are correlated to a lesser extent with ^+ 
at 38 arcmin. The correlation matrix is shown in Fig. 1181 
and the corresponding covariance matrix is tabulated in the 
Appendix in Table E2l 



6.3 Parameter Fits 

We perform all of our fits to a standard five-parameter 
ACDM mode0. For the initial likelihood analysis, we fix Ws, 
0,bh ^ , ^mh^ , and wp at t heir fiducial best-fit WMAP7 val- 



ues (|Komatsu et al.l 201ll ). and vary erg. The upper panel of 



The optical depth to rcionization t is a sixth parameter im- 
plicitly included in the WMAP7 chains, but with no effect on the 
tensing shear correlation function. 



Fig. [19] shows the likelihood of as with all other parameters 
fixed, with a value at the peak and 68 per cent confidence 
interval of O-GSGto lll- 

For a survey of this size and depth, 
the constraints are comparable to the statistically achievable 
confidence limits. 

We also perform a likelihood analysis fixing three pa- 
rameters, and varying f2,„ and erg simultaneously, as these 
two parameters are much more sensitive to the measured 
cosmic shear signal than the others. The resulting two- 
dimensional constraints are shown in the bottom panel of 
Fig. 1191 Our 68 per cent confidence limits on the degen- 
erate product as ( 0^264 )°^^ ^-^^-o.il ^r the Coyote 
Universe prediction code (see Fig. [191 solid red line), and 

as ( 0^264 )°^^ ~ 0-6'''io'i5 foi' the Halofit prediction code (see 
Fig. 1191 dashed red line) . 

We show the effects of removing each systematic error 
correction. Fig. \W\ also shows, for both the one- and two- 
dimensional analyses, the impact of systematic error correc- 
tions. The combined effects of these uncertainties are clearly 
substantially smaller than the statistical error on the ampli- 
tude of the shear signal. 

Finally, we adopt the WMAP7 likelihoods as priors, and 
evaluate our likelihood at each link in the WMAP7 Markov 
chain. For each chain element, we assign a weight equal to 
our likelihood function evaluated at the parameter vector 
for that chain element. For each of the parameter constraint 
plots shown here, we first assign each Markov Chain Monte 
Carlo (MCMC) chain element to a point on a regular grid in 
the parameter space; the value of the marginalised likelihood 
at each grid-point, Hij is then the sum of our likelihood 
weights over the MCMC chain elements at the (i,j) grid- 
point, 



(33) 



where the indicator function Ik{i,j) is equal to unity when 
the grid-point in parameter space is nearest the feth 

chain element, and zero otherwise. The likelihood Lt for 
each chain element is evaluated in the usual way as: 



Lk = exp 



dk 



(34) 



Here C is the full covariance matrix for the measurement, 
incorporating both the statistical and systematic uncertain- 
ties, and the normalization is arbitrary. The data vector dk is 
the extended COSEBI vector described above; where shown, 
the WMAP7 priors are simply this sum with Lk = 1 for each 
point. 

We estimate the detection significance for the final 
signal, the difference 2A log L between the highest- 
likelihood Markov Chain element for both the ACDM and 
uiCDM models and the likelihood evaluated with no sig- 
nal. The Icr detection significances for these two models are 
2.64 and 2.88, respectively. This is not the significance of 
the detection of cosmic shear (as in Sec. 15. 2p . but rather a 
measurement of the likelihood of these two models given the 
combination of WMAP7 priors with this experiment. 

In Fig. 1201 we show marginalized posterior likelihoods 
in the case of fixed ACDM (i.e., w = -1) for Qth'^, 
Us , and as . The results with a free equation of state of dark 
energy (wCDM) are in Fig. 1211 Our measurement provides 
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Figure 19. The effect of systematic errors in tlie 1-D likeliliood of 
(T8 (upper panel) and in tfie 2-D constraints (68 per cent likeliliood 
contours only) in the erg — plane (lower panel). The solid 
curve shows our final analysis, while the other curves show results 
including subsets of the systematic errors. The dot-dashed curve 
labeled "no systematics" shows only the statistical errors, without 
any systematic error corrections either to the theory or to the 
covariance matrix. The dot indicates the WMAP7 best-fitting 
values. 

some additional constraints beyond those from WMAP7 on 
these parameters. In particular, the low amplitude of the 
measured shear signal rules out some of the previously al- 
lowed volume of Q,mh^ and an WMAP7 constraints. 



7 CONCLUSIONS 

Using coadded imaging constructed from SDSS Stripe 82 
data, we constructed a weak lensing catalogue of 1 328 885 
galaxies covering 168 square degrees (Paper I), and showed 
that the additive shear systematics arising from the PSF 
are negligible compared to the cosmic shear signal. In this 
paper, we carried out a cosmic shear measurement that re- 
sulted in a 20 per cent constraint on erg (with all other cos- 
mological parameters fixed). This adds constraining power 
beyond that from WMAP7, and serves as an important 
independent data point on the amplitude of the matter 
power spectrum at late times. In particular, the primary 



CMB anisotropies presently provide only a modest con- 
straint on flrnh^ , and (due to the effect of matter den- 
sity on the growth of structure) there is then an elon- 
gated allowed region in the {Q,mh^,o-s) plane; see Fig. 1201 
The WMAP7-allowed region is ideally oriented for lens- 
ing to play a role: the lensing signal at the high-Qmh^ , 
high-(T8 end of the ellipse leads to a much higher lensing 
signal than low Q,mh? , low (Tg. The low amplitude of cos- 
mic shear observed in this paper eliminates the high-ilm/i^, 
high-erg solutions, and leads to a WMAP7-1-SDSS lensing 
solution of erg = 0.784tj;:°^8(lcr)tE;;°^5(2cr) and n^/i^ = 
0.1303tE!:|5^4g(lcr)j:[];°°^^(2(j); the 2cr error ranges are respec- 
tively 14 and 17 per cent smaller than for WMAP7 alone. 

We have also carefully evaluated other sources of un- 
certainty such as the source redshift distribution, intrinsic 
alignments, and shear calibration, to ensure that our mea- 
surement is dominated by statistical errors rather than sys- 
tematic errors. This achievement is important when consid- 
ering that (i) the SDSS data were never designed with this 
application in mind, and indeed includes several features 
(e.g. the minimal amount of cross-scan dithering) that cause 
significant difficulty, and (ii) with the multitude of upcom- 
ing multi-exposure lensing surveys in the next few years, it 
is important to cultivate new data analysis techniques (such 
as the one used here) that are capable of producing homo- 
geneous data with tight control over PSF anisotropies. As 
a quantitative measure of the extent of PSF correction pos- 
sible with SDSS data, we take the RMS residual spurious 
shear at a particular scale estimated from the star-galaxy 
correlations, 

7rms.oq(6') = • (35) 

From Fig. |6] we see that this is ~ 2 x 10"'' at the small- 
est scales (1-6 arcmin), is < 10"'^ at scales 6 > 0.1 degree, 
and drops to 3.7 x 10"** in the final bin (1.2-2.0 degrees)I!3 
There is almost no difference between the -I— I- and x x sig- 
nals, suggesting that the spurious additive elhpticity signal 
contains similar amounts of E- and B-modeqfj; something 
similar was seen in the SDSS single-epoch data via run-by- 
run comparis ons of elhpticity meas urements on the same 
galaxies ( Ma ndelbaum et all l2006bl . Fig. 8) . This is good 
news for the use of the B-mode as a diagnostic of PSF sys- 
tematics, although an understanding of the generality of this 
pattern remains elusive. 

A major lesson learned from this project is the impor- 
tance of masking bias, in which the intrinsic orientation of a 
galaxy affects whether it falls within the survey mask. This 
is likely the main reason why we had to implement the (ei) 
projection. While we have clearly not exhausted the range of 
options for removing this bias at the catalogue level, future 
surveys should be designed to produce more uniform data 
quality via an appropriate dithering strategy and suppress 
the masking bias at the earliest stages of the analysis. 

Our major limitation in the end was the source number 
density, which was driven by the fact that our PSF-matching 
procedure was limited by the worst seeing in the images that 



" We used TZpsi = 0.9 and Ti = 1.776, as described in the text. 
21 Recall that ^++ (6») - ^ x x (9) and Pe {£) - Pb W are J4 Hankel 
transforms of each other. 
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Figure 20. The cosmological parameter constraints using tlie extended COSEBI data vector, fixing tlie dark energy equation of state w 
at ~1, but allowing all other parameters to vary. Off-diagonal panels show joint two-dimensional constraints after marginalization over 
all the other parameters, which are shown. For these, the red contours show the WMAP7 priors containing 68.5 and 95.4 per cent of 
the posterior probability. The black contours are the same but for WMAP7-(-SDSS lensing. Diagonal panels show the fully-marginalized 
one-dimensional posterior distribution for each parameter; for these panels, the red (dashed) contours show the marginalized WMAP7 
constraints. 



we use, and therefore we had to ehminate the images with 
seeing worse than the median. This means that the coadds 
were not as deep as they could have been, and the final ef- 
fective seeing was 1.31 arcsec (full- width half maximum). In 
principle this will be an obstacle to applying this technique 
in the future, but in fact, that statement depends on context. 
For example, for a survey such as HSC or LSST where we 
expect typically ~ 0.7 arcsec seeing, and with plans to pref- 
erentially use the best-seeing nights for r and i-band imaging 
that will be used for shape measurement, it is conceivable 
that nearly all images intended for lensing will have seeing 
in the 0.6-0.8 arcsec range. In that context, a PSF-matched 
coadd that has the rounding kernel applied may actually 
not result in much loss of information about the shapes of 
most useful galaxies, and will have the advantage of the re- 
moval of PSF anisotropies. Moreover, even for surveys for 
which the loss of information that results from this method 
may not be suitable for the final cosmological analysis, this 
method may still serve as a useful diagnostic of the additive 
PSF systematics. 
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Figure 21. The cosmological parameter constraints using tiie extended COSEBI data vector, varying all five parameters. Off-diagonal 
panels show joint two-dimensional constraints after marginalization over all the other parameters, which are shown. For these, the red 
contours show the WMAP7 priors containing 68.5 and 95.4 per cent of the posterior probability. The black contours are the same but 
for WMAP7-I-SDSS lensing. Diagonal panels show the fully-marginalized one-dimensional posterior distribution for each parameter; for 
these panels, the red (dashed) contours show the marginalized WMAP7 constraints. 
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Table Al. Our data vector. The first five elements are COSEBI 
mode amplitudes; the final is the correlation function averaged in 
the range 29.2296 ^ 6» s: 44.9730. 



4.89797E-10 
1.28335E-09 
1.25136E-09 
1.45616E-09 
8.92333E-10 
1.46457E-05 



Table A2. The covariance matrix for the data vector shown in 
table IAT] 



Data vector index Data vector index Covariance 









3.37161E-20 





1 


4.67637E-20 





2 


4.00484E-20 





3 


2.49916E-20 





4 


9.84257E-21 





5 


3.01770E-17 


1 


1 


1.06383E-19 


1 


2 


1.19226E-19 


1 


3 


8.39508E-20 


1 


4 


3.86519E-20 


1 


5 


1.82344E-16 


2 


2 


1.99923E-19 


2 


3 


1.87469E-19 


2 


4 


1.12196E-19 


2 


5 


5.07790E-16 


3 


3 


2.56568E-19 


3 


4 


2.13363E-19 


3 


5 


8.02118E-16 


4 


4 


2.67774E-19 


4 


5 


5.67797E-16 


5 


5 


3.68112E-11 



APPENDIX A: THE DATA VECTOR AND 
COVARIANCE MATRIX. 

Here we reprint the data vector and covariance matrix used 
in this measurement. The code used to project the correla- 
tion function onto the COSEBI basis functions is available 
from the authors upon request. 
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